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A Particle number conserving shell-correction method 
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The shell correction method is revisited. Contrary to the traditional Strutinsky method, the 
shell energy is evaluated by an averaging over the number of particles and not over the single- 
particle energies, which is more consistent with the definition of the macroscopic energy. In addition, 
the smooth background is subtracted before averaging the sum of single-particle energies, which 
significantly improves the plateau condition and allows to apply the method also for nuclei close 
to the proton or neutron drip lines. A significant difference between the shell correction energy 
obtained with the traditional and the new method is found in particular for highly degenerated 
single-particle spectra (as i.e. in magic nuclei) while for deformed nuclei (where the degeneracy is 
lifted to a large extent) both estimates are close, except in the region of super or hyper-deformed 
states. 

PACS numbers: 21.10.Dr, 21.10.Ma, 21.60.Cs, 21.60.Jz, 25.85.Ca 



I. INTRODUCTION 

The macroscopic-microscopic method of evaluating 
the potential energy surfaces and binding energies of 
nuclei was proposed in the papers of Strutinsky [H 
and Myers and Swi§,tecki 0. Despite the tremendous 
progress of selfconsistent models to nuclear structure the 
macroscopic-microscopic method remains one of the most 
important tools. In such an approach the microscopic 
energy corrections are added to the macroscopic part of 
the nuclear binding energy described by the liquid drop 
model or other macroscopic methods. The microscopic 
part consists of shell and pairing energies. The prescrip- 
tion for the evaluation of the shell energy by smoothing 
the single-particle energy spectra was first given in Ref. 
and than improved in Refs. nil. This Strutinsky 
method of averaging over single-particle energies is still 
widely used up to now, in spite of its known problems 
which appear for mean-field potentials of finite depth as 
well as for nuclei close to the proton or neutron drip lines. 

Already in the 70' (see Refs. |l|E|l|l|l[l3 and re- 
lated papers) Strutinsky and Ivanyuk made an attempt 
to replace the original Strutinsky method of evaluating 
the smooth energy component by an averaging in the 
space of particle numbers (A/'-space) that should be more 
consistent with the macroscopic part of the binding en- 
ergy which is usually evaluated in a liquid-drop type ap- 
proach. The parameters of such macroscopic models are 
usually obtained by a least-square fit to nuclear masses 
which corresponds to an averaging in the A/'-space (e.g. 
in Ref. ^3)- Refs. 0,13 the smooth component of the 
total single-particle energy was approximated by a poly- 
nomial in A/'-space with coefficients that were determined 
by a least-square fit. It was shown in Refs. [1,^3 that the 
shell correction energies obtained by these two types of 
averaging procedures are not the same. Significant differ- 
ences appear for highly degenerated single-particle spec- 
tra, as e.g. in spherical nuclei. The method by Ivanyuk 
and Strutinsky of finding the smooth energy developed 
in Refs. 0, S IE has reached sufficient accuracy to 
be used in practical calculations. It was, however, never 



widely used, probably because of its complexity. 

Another way of separating out the smooth part of the 
sum of single-particle energies can be found in Ref. [T^ . 
where the liquid-drop type asymptotic expansion of the 
total single-particle energy in powers A^^^ was used. Un- 
fortunately this method of evaluating the average energy 
was not precise enough to be used in practice. 

In the present paper a different method of evaluating 
the shell energy is proposed. The smooth component 
of the total single-particle energy is obtained by folding 
the sum of single-particle energies in the A/'-space with a 
modified Gauss function as described in the Appendix. 
In addition, an average energy background as obtained 
by the harmonic oscillator energy sum rule (see section 
III Al below) is subtracted before performing the folding, 
which significantly increases the precision of the method. 
Our new prescription for the shell correction energy gives 
results close to those obtained in the Ivanyuk and Struti- 
nsky approach of Refs. 8, 10] and is extremely simple to 
use. 

One should also mention that the shell energy evalu- 
ated with the present model conserves exactly the given 
number of particles, and not only on the average, as was 
the case in the traditional Strutinsky method. 



II. THEORETICAL MODEL 

In the macroscopic-microscopic method of evaluating 
potential energy one decomposes the nuclear binding en- 
ergy into three parts 

E{Z,A- def) = E,^,,{Z,A- dei) + E,^,n{Z , A- def) 
+EpMZ,A- def) , 

(1) 

where Z and A are the charge and mass numbers re- 
spectively. The macroscopic part, i?mac, depends on the 
deformation of nucleus and is usually evaluated in the liq- 
uid drop or some other more sophisticated model. The 
microscopic part of the energy consists of the shell and 
the pairing energies. The pairing energy E'pair is usually 
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evaluated in the (projected or not) BCS formalism (see 
e.g. 01 or 0), while the shell energy E'sheii is the sum of 
the proton and neutron contributions 

E^iMZ, A; def ) = E^^^^^^iZ; dcf) + E^y^^^.iA - Z; def) . 

(2) 

The shell correction energy of one kind of particles is 
equal to the difference 



(3) 



(10) 



Eq. © leads to the known expression for the average 
density of the harmonic oscillator single-particle levels 



de (S^o)^ 



(3AA) 



2/3 



(11) 



The sum E of single-particle energies of all occupied lev- 
els is 



where Af is the number of particle in the system and E is 
the smooth part of the total single-particle energy, where 
smooth means slowly varying with the particle number 
Af. In the following two different methods of evaluating 
of this smooth part will be presented. 



A. Harmonic oscillator energy sum rule 
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E = Y, en deg„ = ^c^^ ^ (n + -) (n + 1) (n + 2) (12) 

n=0 n=0 

and can be approximated by the integral 

AT 

E w / e(AA') dN' . (13) 



The eigenenergies of the spherical harmonic oscillator 



e„ = (n + -) hojo 



are strongly degenerated 
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deg„ = -(n+l)(n + 2) x 2 



(4) 



(5) 



Here n = 0, 1, 2, ... is the main quantum number and loq 
is the harmonic oscillator frequency. The factor 2 in the 
above equation is due to the two possible orientation of 
the spin. 

According to Ref. jl3| the degeneracy of the main har- 
monic oscillator shell can be approximated by 



o\ 2 / \ 2 



(6) 



The total number of particles Af occupying all shells up 
to n = N is 

AA(N) = ^deg„ = -(N-f l)(N + 2)(N-f 3) . (7) 

n=0 

It is easy to show ^3 that for large N values the following 
approximation holds: 



(8) 



The last equation can serve as the average relation be- 
tween the single-particle energy e and the number of 
particles which occupy the levels with energy smaller or 
equal to e 



(9) 



Inserting here Eq. H1(J|I one obtains the following energy 
sum rule: 



(14) 



The sum of energies of nucleons which occupy the har- 
monic oscillator levels is thus proportional to the 4/3 
power of the total number of particles in the system. 
A more accurate estimate than the above one was made 
in Ref. m : 
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(15) 



The term proportional to Af^^^ is important in the light 
systems but in the heavier nuclei it can be neglected as 
much smaller than the leading Af^^^ term. 

In the top l.h.s. part of Fig. ^ the sum E of single- 
particle energies (solid line) and its approximation E 
(dashed line) by Eq. 1)14(1 are shown as function of the 
number of particles Af. The deviation AE between both 
lines is hardly visible on this scale, so we present it sep- 
arately in the top r.h.s. part of Fig.Q] The coefficient in 
front of the term Af^^^ was obtained by a least square fit 
and turns out to be very close to the value of the approx- 
imate expression (I14|l which is exact in the limit N — > oo. 
A strong shell structure corresponding to the harmonic 
oscillator magic numbers: 7V„= 2, 8, 20, 40, 70, 112, 168, 
240, 330, 440, .... is observed. Using Eq. © one can ob- 
tain the average distance between the harmonic oscillator 
major shells as function of the particle number 



(16) 



The deviation AE of the energy sum from its average 
behavior as function of Af^^^ is presented in the bottom 
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FIG. 1: Sum of the single-particle energies {E, solid line) and its approximation {E, dashed line) by Eq. I|14|l (top l.h.s. plot) 
as well as their difference {AE — E — E) in function of the number of particles Af (top r.h.s. plot), 7V^^^ (bottom l.h.s. plot) 
and the single particle energies (e^) (bottom r.h.s. plot). 



l.h.s. part of Fig. ^ It is seen that the distance be- 
tween closed shells is nearly constant and roughly equal 
to A(7V^/^) w 0.7 which is the estimate, Eq. It is 

worth noticing that the same data plotted as function of 
the single-particle energies e shows a structure (bottom 
r.h.s. of Fig. ^) which seems hard to interpret at first 
sight. Obviously the shell structure of the harmonic os- 
cillator is more visible when one plots AE as function of 

The relation (|14|l was obtained assuming that the 
single-particle energies are measured with respect to en- 
ergy zero. Assuming that the minimum of the harmonic 
oscillator potential corresponds to Vq (i.e. et ^ Ci + Vq) 
one can get the more general relation 



^ = (^E ~ + Vo A/- . (17) 

We have verified (numerically) that the above har- 
monic oscillator energy sum rule is universal and not 
only fuUfiled by the spectra of the modified harmonic 
oscillator (Nilsson potential) or other finite depth model 



mean-field potentials (e.g. Saxon- Woods) but also by 
the single-particle spectra obtained selfconsistently for 
the Hamiltonians associated with the Gogny or Skyrme 
effective forces. 

A typical deviation of the sum of the single-particle en- 
ergies (with respect to the bottom of the effective mean- 
field potential) from the estimate (|17|l is of the order of 
a few promilles for heavier nuclei. In Fig. |21 the sum of 
the single-particle energies (l.h.s. column) and its de- 
viation (AE) (r.h.s. column) from the average trend, 
Eq. (|17|l . is plotted as function of Z'^^^ for protons (top 
row) and N^^^ for neutrons (bottom row). The single 
particle energies of spherical ^'^^Pb were obtained self- 
consistently using the Hartree-Fock approximation to the 
Gogny Hamiltonian with the DIS force T^. The param- 
eters a and Vq of Eq. H17|) given in Fig.[21are obtained by 
a least square fit. The arrows point the Fermi-level po- 
sitions and the dotted vertical lines mark the end of the 
bound state spectrum. A very pronounced shell struc- 
ture of the proton and neutron spectra is visible in the 
r.h.s. plots. 
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FIG. 2: Sum of the single-particle energies {E, solid line) obtained selfconsistently with the Gogny DIS force for ^°^Pb and its 
approximation (E, dashed line) by Eq. H17|l as well as the deviation from this average trend (r.h.s. plots) as function of the 
number of protons (top) or neutrons (bottom). Arrows indicate the position of the Fermi energy and the vertical lines mark 
the end of the bound spectrum. 



B. Average of the sum of single- particle energies 

Let us define a discrete sample of data, Sn, as the 
difference between the sum of ttie lowest available single- 
particle energies of the n fermion system and the corre- 
sponding background energy, E{ti), obtained using the 
harmonic oscillator sum rule, Eq. H17|) 

n n 

Sr, = Y,e,-E{n)=Y,e,^an^'^ -V^n . (18) 
1=1 1=1 

The parameters a and Vq are determined by minimizing 
the square deviation between the single-particle energy 
sum and E 

n=l 

where Mmax can be chosen as the maximal number of 
nucleons which can be put on the given single-particle 
energy spectrum. 



Using the Gauss-Hermite folding procedure described 
in details in Appendix A one can evaluate the average 
value of Sn corresponding to M nucleons 




where /g is the 6*'* order polynomial given by Eq. ljA26|l . 
The folding is performed not directly in the particle num- 
ber n but in its cubic root since the distance between the 
major harmonic oscillator shells is constant in 'n}/'^ and 
approximately equal 0.7 as we have shown above. The 
factor 3n^/'^ in the denominator of Eq. (|20|l is the direct 
consequence of the transformation n n^/'^, while the 
factor 2 in the numerator is due to the spin degeneracy 
of the single-particle levels. 

The smoothed energy of an even or odd A/" system is 
then 

E{M) ^Sm + oM^'^ + VqM , (21) 
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where we have restored the background energy E(Af), 
Eq. H17|l . which has been subtracted from the single- 
particle energy sum in Eq. H18|l . Subtracting E{n) in 
(|18|l increases significantly the accuracy of evaluating the 
smoothed part E of the energy since the deviations S'„ 
is 2 to 3 orders of magnitude smaller than the value of 
E{M). The smoothed energy obtained in this way is less 
sensitive to the energy cut-off of the single-particle spec- 
trum, which is important for evaluating the shell energy 
of nuclei close to the proton or neutron drip lines. 

In Fig. O the sum of the single-particle energies (solid 
line) is compared with the new particle-number smoothed 
energy (dashed line) and the old Strutinsky energy (dot- 
ted line). The single particle spectrum is the one evalu- 
ated for the spherical Saxon- Woods mean-field potential 
for the ^'^^Pb nucleus with the parameters taken from 
Ref. m. The background energy E is subtracted from 
the all three energies presented separately for protons 
(l.h.s. plot) and neutrons (r.h.s. plot). One notices that 
the Strutinsky energy is always smaller than the present 
estimate of the smoothed energy. The difference between 
both estimates grows with the number of particles and for 
208p^ (arrows in Fig.|21) is of the order of a few MeV. This 
result is similar to the one obtained with the Ivanyuk and 
Strutinsky method ^3 j where the smoothed part of 
the sum of single-particle energies was approximated by 
a local polynomial in the A/'-space. 



C. Strutinsky smoothed energy 

It is worthwhile to remind here the original Strutin- 
sky method of evaluating of the smooth energy ^l, jl, ^ 
in order to better understand the difference between 
both approaches. The Strutinsky 's way of evaluating the 
smooth energy consists of two steps. Firstly one finds 
the smoothed single-particle level density 5(e) and deter- 
mines the corresponding average position A of the Fermi 
level, assuming the average particle number conservation. 
Then in the next step one evaluates the smoothed energy 
by integrating the product of the single-particle energy 
and smooth level density. It means that in this method 
the number of particles is conserved only on the aver- 
age and the Strutinsky smoothed energy does not corre- 
spond exactly to the averaged sum of the occupied single- 
particle energies. 

In the Strutinsky shell correction method one evaluates 
the smooth single-particle level density g{e) by folding 
the discrete spectrum of eigenstates 



5(e) = ^<5(e-ei) 



(22) 



with a smoothing hmction jVi(e, e ) of the n*'* order which 
is given by Eq. (|A11|) . The smooth single-particle level 
density 5(e) is then given by 



5(e) = y^jn { - — - 



(23) 



Taking the 6*^ order (so called "curvature corrections^) 
polynomial into account (see Eq. (|A26|l ) the smoothing 
function has the following form 



,,2,35 35 

" ( 1 
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(24) 



The smearing parameter 75 in Eqs. (|23l I24|l is the width 
of the Gauss folding function and should be of the order 
of the energy distance between major shells (i.e. hcoo) in 
order to wash out the shell structure. 

According to Strutinsky Q the average of the sum of 
the energies of the occupied single particle levels (i?sti ) 
is given by the integral 



£^str = J 2eg{e)de , 



(25) 



where A is the position of the Fermi energy in the system 
with the washed out shell structure and is fixed by the 
particle number condition 



2 g{e) de 



(26) 



Here the average number of particles Af = Z ior pro- 
tons or M = N for neutrons. The factor 2 in the above 
two equations is due to the spin degeneracy of the single 
particle levels. One solves Eq. H26(l for A by iterations. 

The Strutinsky energy Estr l|25(l is not equal to the av- 
erage of the sum of single-particles energies E, Eg. I21|l . 
but corresponds to the energy of a system which con- 
serves the number of particles only on the average (and 
not exactly as in Eq. H21|l ). A comparison of the result- 
ing smoothed energies obtained in both methods will be 
presented below. 



III. COMPARISON OF THE BOTH ESTIMATES 
OF THE SMOOTHED ENERGY. 

A significant difference between the new estimate of the 
smooth energy E given by Eq. (|21|l and the Strutinsky 
energy, E'str, Eq. (|25|) . is demonstrated in Fig. 13 It can 
be easy explained as follows : 

The sum of single-particle energies can be roughly ap- 
proximated by Eq. lfn|) as shown in Sec. Ill Al 




(27) 



where the parameter a is proportional to the distance 
hujQ between major shells and b to the effective depth of 
the mean-field potential. Let us assume, just as a matter 
of discussing our method, that this average trend repre- 
sents the true energy sum and that we are dealing with 



6 




FIG. 3: Sum of the single-particle energies E (solid lines) obtained for the Saxon- Woods potential of Pb and its smooth part 
obtained with Eq. I21II (fat dashed lines) as well as within the traditional Strutinsky method (thin dotted lines). From all the 
curves is subtracted the background energy evaluated as in Eq. 11711 . The data for protons and for neutrons are presented in 
the left and right parts respectively. Arrows indicate the positions of the Fermi energies and the vertical dotted lines indicate 
the end of the bound spectra. 



the degenerate spectrum. Let also A4 be the number of 
particles which can be placed on the single-particle levels 
which are below the k-th degenerate level. Note, that 
the numbers A4 (with fc = 1, 2, ...) are simply the magic 
numbers in case of spherical nuclei. 

The Strutinsky prescription for the smoothed energy 
corresponding to 

^f = ^i^fk+i - ^fk) (28) 

particles (i.e. half filled shell) can be written as 

Estr = aM^'^ + bM , (29) 

while the average over the interval [WfejAfc+i] of the en- 
ergy, Eq. (|27|l is given by the integral 



E 



1 

A 



A/'+A/2 



{an^^^ + bn) dn , 



(30) 



AA-A/2 



where A = Afk+i — Mk is the degeneracy of the corre- 
sponding single-particle level. In the approximation H27|l 
the diS'erence between the average of the single-particle 
energy sum and the Strutinsky energy is 



AE = E-Ei 



Str 



1 



54 



(31) 



where a « 3'*^'^/4 fiwo- This approximate expression in- 
dicates that the diS'erence between the new and the old 
Strutinsky energy is negligible (of the order of 0.01 MeV 
for heavier nuclei) when the degeneracy is A = 2 which is 
the case for deformed nuclei, while it grows significantly 
(up to a couple of MeVs) when the degeneracy is impor- 
tant, as e.g. in spherical nuclei. 




20 40 60 80 100 120 140 160 180 

N 

FIG. 4: Strutinsky shell energies for a spherically symmetric 
harmonic oscillator potential obtained with the new (Eq. 1211 
solid line) and the old f Eg. 1251 dashed line) method (see text) 
as well as the difference between them (dotted line) as func- 
tion of number of particles A/". The estimates for the differ- 
ences AE, given by Eq. ^ are marked as the full circles. 



In Fig. 0] we show the shell energy EshcW for the spher- 
ical harmonic oscillator single-particle levels in function 
of the nucleon number (Af) . The solid line corresponds to 
the new approach described in Sec. lIIBl while the dashed 
one to the old Strutinsky method (see Sec. Ill Cp . The 
difference between the both estimates and its approxima- 
tion with Eq. (|31|l are shown by the dotted line and the 
full circles respectively. 

The proton and neutron shell energies -EshcU for the 
nuclei ^^^Pb and ^'^^Th obtained with a deformed Saxon- 
Woods potential are plotted in Fig. as function of 



7 




FIG. 5: Proton (left) and neutrons (right) Strutinsky shell energies for the Saxon- Woods single-particle levels of ^°^Pb (top) and 
232 rpj^ (bottom) obtained with the new (Eg. 1211 solid line) and the old (Eg. 1251 dashed-line) method as well as their difference 
(dotted line) as functions of the elongation parameter c 3- 



the elongation parameter c of Ref. j^]. The sohd line 
corresponds to the new approach described in Sec. Ill Bl 
while the dashed one to the old Strutinsky method (see 
Sec. Ill C|l . The difference between the both estimates is 
shown by the dotted line. The parameters of the Saxon- 
Woods potential are taken from Ref. (0|). It is seen 
that the difference between the shell energies evaluated 
with the new particle number conserving method and the 
traditional Strutinsky approach becomes negligible with 
growing nuclear deformation, i.e. when the degeneracy 
of single particle levels is lifted, a result which is in line 
with the prediction of the approximate expression (|31|l . 

In both (new and old) Strutinsky shell correction meth- 
ods it is very important to choose the appropriate value 
of the smearing parameter (7 or 75). This is usually done 
by fixing its value such that the obtained shell energy is, 
over a certain range of that parameter independent of its 
specific value, a constraint known as the ^^plateau con- 
dition^^ of the Strutinsky method. Typical examples of 
such plateaus are shown in Fig. ^ where the shell ener- 
gies for ^"^Pb at a deformation of c = 1.2 obtained with 
the appropriately chosen Saxon- Woods potential J^] are 



drawn. It is seen in Fig. that our new approach leads 
to a much more pronounced plateau as compared to the 
old method. Note the different values of the smoothing 
parameter (7 « 0.78 for both kinds of particles versus 
7s ~ 1.20?ia;o for protons and 7s ~ l.OdhuJo for neu- 
trons). 

Looking at Fig. |S1 one could have the impression that 
a significant difference between the shell energies ob- 
tained by the averaging in the particle-number or single- 
particle energy spaces appears around spherical shapes 
and vanishes with growing deformation. On the other 
hand, Eq. H31|l predicts that a significant difference be- 
tween both types of averaging procedures should appear 
whenever a large degeneracy of the single-particle levels 
is present. A good example of the spectrum which be- 
comes strongly degenerate at some deformation points 
are the eigenenergies of the anisotropic harmonic oscilla- 
tor. The degeneracy clearly appears at those ellipsoidal 
deformation (e, Q) points where the ratio of the axes 
is equal to the ratio of small integers as can be seen in 
the left part of Fig. [7| where such a spectrum is plotted. 
The harmonic oscillator single-particle levels presented 
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FIG. 6: Dependence of proton (solid line) and neutron (dashed line) Strutinsky shell energies as function of the smearing 
parameter 7 (or 7s) obtained with the new (Eq. I21II . left) and the old (Eq. I25II . right) estimates for the smooth energy for a 
nucleus ^°*Pb described by a appropriate Saxon- Woods potential at a deformation of c = 1.2. 



10 
9 
8 
7 

3" 6 

4 
3 
2 
1 



_l I I I I I l_ 



-0.5 




—I I I l_ 



0.5 



2 


1 1 1 1 1 1 1 1 1 
Harmonic oscillator 


' 1 ' 
1:1 


1 1 1 1 1 1 1 1 1 
W= 100 


1.5 








3 


1:2 






LJJ 






2:1 


LU 

0.5 


-2:5 / \ 




l\ 3:1 
3:2 , 





" 1 , , , , 1 , , , 


, 1 , 


, , , 1 , , , , 1 ~ 




-1 -0.5 




E 


0.5 1 



FIG. 7: The eigenfunction of the deformed harmonic oscillator potential (l.h.s. plot) and the difference between the new(E, 
Eq. I2II1 ) and the old {Estr, Eq. 12511 ) estimates of the smooth energy as function of the axial quadrupole deformation e Q. 



there were used to evaluate the sheh energy of a system 
composed of J\f=100 particles. The difference between 
the new and the old estimates {E and E'str) is plotted 
as function of the quadrupole deformation parameter e 
[3- A large difference between both approaches appears 
where the degeneracy of levels grows. A similar effect was 
already observed in Ref. 0- This result clearly shows 
that using the new approach one can expect some modi- 
fications in the potential energy surface not only around 
spherical shapes but also at deformations which corre- 
spond to the super- or hyper- deformed isomers. 



IV. SUMMARY AND CONCLUSIONS 

A new method of evaluating the smooth part of the 
total single-particle energy is proposed. The folding of 



the sums of single-particle energies is performed in the 
particle-number space rather than in the one of single- 
particle energies as done in the old Strutinsky method. 
The averaging in the A/'-space is consistent with the def- 
inition of the macroscopic energy component which rep- 
resents the average behavior in Z and A of the nuclear 
binding energy. 

One has also to notice that the integral over Af of the 
shell energy evaluated with the new prescription is close 
to zero for sufhciently large number of particles, while 
this is not the case for the traditional Strutinsky shell 
correction which grows systematically with J\f. This indi- 
cates that the original Strutinsky prescription generates 
shell energy which do not fluctuate around zero when the 
number of particles is increased which modifies system- 
atically the macroscopic part of the energy. In addition 
this modification depends on the shape of nucleus, so 
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that this deficiency of the traditional approach cannot 
be corrected by an adjustment of the parameters of the 
macroscopic energy. 

The new estimate differs significantly from the Struti- 
nsky smoothed energy when a large degeneracy of the 
single-particle levels is present as this is the case in spher- 
ical and nearly-spherical nuclei or in some shape iso- 
mers. In such nuclei the shell energy is shifted down by 
a few MeV with respect the old predictions while its am- 
plitude is almost unchanged as function of the particle 
number. This means that the macroscopic-microscopic 
method with the new estimate of the shell energy will, 
when leaving the parameters of the macroscopic and mi- 
croscopic parts untouched, predict the magic and quasi- 
magic nuclei as well as some shape isomers more bound 
than this was predicted by the calculations done with the 
old Strutinsky method. Also the deformation energies of 
non-magic nuclei will be smaller and one could obtain dif- 
ferent equilibrium shapes (i.e. ground state quadrupole 
moments). The fission barrier for spherical nuclei will 
be also significantly increased and, as a consequence, the 
spontaneous fission of such nuclei (e.g. some super-heavy 
isotopes) will be less probable. In addition the Q-value 
for an a-decay will be modified when it occurs between 
deformed and spherical isotopes (or vice versa). This 
means that the consequences of the naive use of new 
method could be dramatic. 

I would therefore like to end with the following warn- 
ing: 

Do not use the new prescription for the shell energy (Eqs. 
hl^) - in practical calculations without an appro- 

priate readjusting of the parameters of the mean-field 
potentials (e.g. Saxon-Woods, Nilsson, Yukawa-folded), 
macroscopic models (e.g. liquid drop, finite range droplet 
or Thomas- Fermi) , and the pairing force. 
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APPENDIX A: FOLDING OF DISCRETE DATA 

1. General formulae 

Our aim is to approximate a sample of N ordered 
points {xi,yi\ by a continues smooth function y{x). We 
would like to solve this problem using the Gauss-Hermite 
folding method which was idea the originally proposed by 
V.M. Strutinsky Q and later-on generalized in Ref. [3j. 



Having the width of the folding function comparable with 
the average distance between points Xi one can obtain the 
folded function which is very close to the data points but 
increasing this width one can also wash out the fine struc- 
ture stored in the data. Usually the Strutinsky method 
was used to realize the second scope. The parameter of 
the folding procedure will be determined by requirement 
that the integral of the folded function should be the same 
as the integral evaluated with the sample of {xi, yi\ pairs 
using the trapezium rule. 

Let jn{x,x') be a symmetric function of its arguments 
(i.e. jn{x,x') = jn{x' ,x)) having the following proper- 
ties: 

-f-oo 

j jn{x,x')dx = 1 (Al) 

— OO 

and 

+ 00 

Pkix)^ J Pk{x')jn{x,x')dx' , (A2) 

— oo 

where k < n are even natural numbers and Pk{x) is an 
arbitrary polynomial of order k. In the following, the 
function jn{x,x') will be called the folding function of 
the n"^ order. The last equation , frequently called the 
Strutinsky condition ensures that the folding does not 
change the average behavior of the function Y(^x) which 
is represented by the ensemble of {xi,yi} points. An 
example of such a folding function can be a combination 
of the Gauss function and the Hermite polynomials of the 
argument proportional to jx — cc'l, frequently used in the 
Strutinsky shell correction method 0, 0| . More detailed 
description of such a folding function will be given in the 
next section. 

With each discrete point {xi,yi) one can associate the 
function yi{x) defined by: 

+ C30 

y,{x) = J yi 6{x' - Xi) jn{x, x') dx' , (A3) 

— oo 

where 5{x) is the Dirac 5- function. A straightforward 
calculation gives 

yi{x) = yijn{x,Xi) . (A4) 

Using Eq. HA1|I it is easy to verify that the integral of the 
function yi{x) is 

-\-oo 

J yi{x) dx = yi . (A5) 

— oo 

Let us construct the function y{x) by summing up, with 
weight Wi all functions yi{x) 

N 

y{x) ^^w^yi{x) . (A6) 

i=l 
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The function y{x) is an approximation of y{x) if the 
weights Wi are determined from the assumption that the 
integrals of the unfolded and folded function are (nearly) 
equal: 

N +?° N 

^y{xi)Axi^ / yix)dx = ^Wiy, , (A7) 

j=l _^ i=l 



where Axi is set to: 



Ax^ = - {xi+i - x,^i] 



(A8) 



Eq. (|A7(I implies that a reasonable choice of the weight 
is 

= Axi . (A9) 
Thus the folded function y{x) is given by 

N 

y{x) = ^ 2/i Axi j„(x, Xi) . (AlO) 

i=l 

2. Gauss-Hermite folding function 

Let the folding function jn{x,x') be defined with the 
help of the Gauss function as 



in {x^X ^ 



1 



: exp 



(All) 

where 7 is a parameter and fni. ^ ^ ) is the so called 
corrective polynomial of n*'' order, determined by the 
Strutinsky condition l|A2p . In the following we would 
like to evaluate the coefficients of the corrective polyno- 
mial using some properties of the Hermite polynomials 
which are orthogonal with the weight equal to the Gauss 
function. 

Let us introduce variable u= {x~x')/^ defined in the 
interval (—00, +cx)). The smearing function jn{x, x') and 
the polynomial P„(x) in ljA2p can now be written as 



fn{u) 



(A12) 



Pn{x)=Pn{x-lu) = Pj{u), (A13) 



and 



/'„(x) =F„(x + 7 0) = P„'(0). 



(A14) 



The so far arbitrary polynomial {u) can be written 
down as a series of the Hermite polynomials of order i 



P„'(u) =^a,i7,(7/) 



(A15) 



Now the condition (jA2|) can be written as 

+00 



P,/(0) 



1 



/TT 



Pn'(")e"" fn{u)du (A16) 



and inserting relation l|A15p into (|A16p one obtains 

Y^a,^-j=j e-^ H,{u) fr,{u) du - H,{0)\ = Q . 

(A17) 

On the other hand, the last equation should be fuUfiled 
for arbitrary values of at 7^ 0) what leads to the following 
set of equations 

r e-"' [u) /„ {u) du^H,{0) , (A18) 

where i = 0, 2, ...,n. From the other side the corrective 
function /„(u) can be also decomposed in terms of the 
Hermite polynomials 



(A19) 



fc=i 



Inserting the above relation into Eq. IjAlSp gives 



1 I' + OO 



e~'''H,{u)Hk{u)du . (A20) 



Then using the orthogonality properties of the Hermite 
polynomials 



1 



+00 



e-" H,{u) Hk{u) du = T i\ 5,k , (A21) 



one obtains the coefficients of the corrective polynomial 

a = ^ff»(0) (A22) 
The values of the Hermite polynomials at zero are 

1 for i = , 

Hi{Q) = { 2" (-l)"(2n - 1)!! for i = 2n , 

for i = 2n + 1 , 



so that 



(A23) 



for 1 = 0, 



C^={ (-1)"^^ for* = 2n>0, (A24) 







for i = 2n + 1 



The first few coefficients Ci and the corresponding Her- 
mite polynomials are: 

Co - 1 , Hoiu) = 1 , 
C2--i , H2(") = 4u2-2 , 
C4 = 4-^ , Ha{u) = l&u^ - 48u2 + 12 , 
Ce = -3^4 , Heiu) = 64^^ - 480^^^ + 720^^ - 120 , 

(A25) 
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and the resulting corrective polynomials have the follow- 
ing form 



foiu) 
U{u) 

/6M = i-f«^ + K-^^ 



1 , 

3 _ 
2 
15 
8 



(A26) 



Finally the function y{x) approximated using the 
Gauss-Hermite folding reads: 



In principle the smearing parameter 7 is arbitrary and 
it can be different at each point Xi. But it should be 
related to the distance Axi between subsequent points if 
one would like to approximate the function stored in the 
mesh of {xi,yi} points. Similarly one has to choose 7 of 
the order of the period-length of the fine structure (e.g. 
shell effects) in case when one would like to wash out this 
structure from the function y{x). 



1 ^--v ( X — Xi\ J. f X — Xi\ 

^-^g,,A..exp -(-^) 



7 

(A27) 
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